%run preparation.py
%matplotlib inline
HTML('''<script>
code_show=true;
function code_toggle() {
if (code_show){
$('div.input').hide();
} else {
$('div.input').show();
}
code_show = !code_show
}
$( document ).ready(code_toggle);
</script>
<h3>Actin grinding velocity analysis</h3><br># Prepared by Sam., Prof. Nitta Lab.<br># Click <a href="javascript:code_toggle()">HERE</a> to toggle on/off the raw input code.''')
Dx_tip = np.diff(bmd01['x']); Dy_tip = np.diff(bmd01['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(10*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.around(np.linspace(0,3,30,endpoint=True), decimals=1)
#time = np.around(time, decimals=1)
#print(time)
fig, ax1 = plt.subplots(figsize=(10,6), dpi=500)
ax1.set_xlabel('Time (sec)', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm01, label='Change in binding motor')
#ax1.axhline(float(bm01.mean()), linestyle='--')
ax1.legend(loc='upper left')
ax1.set_ylim(bottom=0)
ax2 = ax1.twinx()
ax2.set_xlabel('Time (sec)', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')
ax2.set_ylim(bottom=0)
plt.title('Binding motors [R = 0.1]', fontdict=font)
plt.grid()
v1 = v.reshape(30,1)
v1 = pd.DataFrame(v1)
tym = pd.DataFrame(time)
bm01_ = bm01.reset_index(drop=True)
nbm01 = pd.concat([tym, bm01_, v1], axis=1)
nbm01.columns=['time','b_m','vel']
plt.show()
fig.savefig('actin_binding_grinding_r0.1.svg',bbox_inches='tight', format='svg')
print(colored('Ave_vel:', 'yellow', attrs=['reverse', 'blink']), \
colored(Av_vel, 'yellow', attrs=['reverse', 'blink']))
print(colored('VSD:', 'yellow', attrs=['reverse', 'blink']), \
colored(vSD, 'yellow', attrs=['reverse', 'blink']))
print(nbm01)
Dx_tip = np.diff(bmd02['x']); Dy_tip = np.diff(bmd02['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(10*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.around(np.linspace(0,3,30,endpoint=True), decimals=1)
v1 = v.reshape(30,1)
v1 = pd.DataFrame(v1)
tym = pd.DataFrame(time)
bm02_ = bm02.reset_index(drop=True)
nbm02 = pd.concat([tym, bm02_, v1], axis=1)
nbm02.columns=['time','b_m','vel']
nbm02 = nbm02.loc[nbm02['b_m'] >= 20] # detached number
nbm02 = nbm02.iloc[:,2] # pick column v1
#print(nbm02)
nbm02 = nbm02[0:5]
nv = nbm02; nAv_vel = nbm02.mean()
nvSD=np.sum(((nv-Av_vel)**2)/(np.size(nv)-1)); nvSD=np.sqrt(nvSD)
fig, ax1 = plt.subplots(figsize=(10,6), dpi=500)
ax1.set_xlabel('Time (sec)', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm02, label='Change in binding motor')
#ax1.axhline(float(bm01.mean()), linestyle='--')
ax1.legend(loc='upper left')
ax1.set_ylim(bottom=0)
ax2 = ax1.twinx()
ax2.set_xlabel('Time (sec)', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(nbm02.mean(), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')
ax2.set_ylim(bottom=0)
plt.title('Binding motors [R = 0.2]', fontdict=font)
plt.grid()
plt.show()
fig.savefig('actin_binding_grinding_r0.2.svg',bbox_inches='tight', format='svg')
#print(colored('Ave_vel:', 'yellow', attrs=['reverse', 'blink']), \
# colored(Av_vel, 'yellow', attrs=['reverse', 'blink']))
#print(colored('VSD:', 'yellow', attrs=['reverse', 'blink']), \
# colored(vSD, 'yellow', attrs=['reverse', 'blink']))
print(colored('Mean velocity, time[0:0.4] = ', 'yellow', attrs=['reverse', 'blink']), \
colored(nbm02.mean(), 'yellow', attrs=['reverse', 'blink']))
print(colored('bVSD:', 'yellow', attrs=['reverse', 'blink']), \
colored(nvSD, 'yellow', attrs=['reverse', 'blink']))
nbm02 = pd.concat([tym, bm02_, v1], axis=1)
nbm02.columns=['time','b_m','vel']
print(nbm02)
Dx_tip = np.diff(bmd03['x']); Dy_tip = np.diff(bmd03['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(10*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.around(np.linspace(0,3,30,endpoint=True), decimals=1)
v1 = v.reshape(30,1)
v1 = pd.DataFrame(v1)
tym = pd.DataFrame(time)
bm03_ = bm03.reset_index(drop=True)
nbm03 = pd.concat([tym, bm03_, v1], axis=1)
nbm03.columns=['time','b_m','vel']
#print(nbm03)
nbm03 = nbm03.loc[nbm03['b_m'] >= 20] # detached number
nbm03 = nbm03.iloc[:,2]
nbm03 = nbm03[0:2]
#print(nbm03)
nv = nbm03; nAv_vel = nbm03.mean()
nvSD=np.sum(((nv-nAv_vel)**2)/(np.size(nv)-1)); nvSD=np.sqrt(nvSD)
fig, ax1 = plt.subplots(figsize=(10,6), dpi=500)
ax1.set_xlabel('Time (sec)', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm03, label='Change in binding motor')
#ax1.axhline(float(bm01.mean()), linestyle='--')
ax1.legend(loc='upper left')
ax1.set_ylim(bottom=0)
ax2 = ax1.twinx()
ax2.set_xlabel('Time (sec)', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(nbm03.mean(), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')
ax2.set_ylim(bottom=0)
plt.title('Binding motors [R = 0.3]', fontdict=font)
plt.grid()
plt.show()
fig.savefig('actin_binding_grinding_r0.3.svg',bbox_inches='tight', format='svg')
#print(colored('Ave_vel:', 'yellow', attrs=['reverse', 'blink']), \
# colored(Av_vel, 'yellow', attrs=['reverse', 'blink']))
#print(colored('VSD:', 'yellow', attrs=['reverse', 'blink']), \
# colored(vSD, 'yellow', attrs=['reverse', 'blink']))
print(colored('Mean velocity, time[0:0.1] = ', 'yellow', attrs=['reverse', 'blink']), \
colored(nbm03.mean(), 'yellow', attrs=['reverse', 'blink']))
print(colored('bVSD:', 'yellow', attrs=['reverse', 'blink']), \
colored(nvSD, 'yellow', attrs=['reverse', 'blink']))
nbm03 = pd.concat([tym, bm03_, v1], axis=1)
nbm03.columns=['time','b_m','vel']
print(nbm03)
Dx_tip = np.diff(bmd04['x']); Dy_tip = np.diff(bmd04['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(10*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.around(np.linspace(0,3,30,endpoint=True), decimals=1)
v1 = v.reshape(30,1)
v1 = pd.DataFrame(v1)
tym = pd.DataFrame(time)
bm04_ = bm04.reset_index(drop=True)
nbm04 = pd.concat([tym, bm04_, v1], axis=1)
nbm04.columns=['time','b_m','vel']
#print(nbm04)
nbm04_1 = nbm04[0:4]
nbm04_2 = nbm04[14:20]
frames04 = [nbm04_1,nbm04_2]
nbm04_3 = pd.concat(frames04)
nbm04 = nbm04_3.iloc[:,2]
nv = nbm04; nAv_vel = nbm04.mean()
nvSD=np.sum(((nv-nAv_vel)**2)/(np.size(nv)-1)); nvSD=np.sqrt(nvSD)
fig, ax1 = plt.subplots(figsize=(10,6), dpi=500)
ax1.set_xlabel('Time (sec)', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm04, label='Change in binding motor')
#ax1.axhline(float(bm01.mean()), linestyle='--')
ax1.legend(loc='upper left')
ax1.set_ylim(bottom=0)
ax2 = ax1.twinx()
ax2.set_xlabel('Time (sec)', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(nbm04.mean(), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')
ax2.set_ylim(bottom=0)
plt.title('Binding motors [R = 0.4]', fontdict=font)
plt.grid()
plt.show()
fig.savefig('actin_binding_grinding_r0.4.svg',bbox_inches='tight', format='svg')
#print(colored('Ave_vel:', 'yellow', attrs=['reverse', 'blink']), \
# colored(Av_vel, 'yellow', attrs=['reverse', 'blink']))
#print(colored('VSD:', 'yellow', attrs=['reverse', 'blink']), \
# colored(vSD, 'yellow', attrs=['reverse', 'blink']))
print(colored('Mean velocity, time[0:0.3, 1.4:2.0] = ', 'yellow', attrs=['reverse', 'blink']), \
colored(nbm04.mean(), 'yellow', attrs=['reverse', 'blink']))
print(colored('bVSD:', 'yellow', attrs=['reverse', 'blink']), \
colored(nvSD, 'yellow', attrs=['reverse', 'blink']))
nbm04 = pd.concat([tym, bm04_, v1], axis=1)
nbm04.columns=['time','b_m','vel']
print(nbm04)
Dx_tip = np.diff(bmd05['x']); Dy_tip = np.diff(bmd05['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(10*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.around(np.linspace(0,3,30,endpoint=True), decimals=1)
v1 = v.reshape(30,1)
v1 = pd.DataFrame(v1)
tym = pd.DataFrame(time)
bm05_ = bm05.reset_index(drop=True)
nbm05 = pd.concat([tym, bm05_, v1], axis=1)
nbm05.columns=['time','b_m','vel']
#print(nbm05)
nbm05 = nbm05[0:4]
nbm05 = nbm05.iloc[:,2]
nv = nbm05; nAv_vel = nbm05.mean()
nvSD=np.sum(((nv-nAv_vel)**2)/(np.size(nv)-1)); nvSD=np.sqrt(nvSD)
fig, ax1 = plt.subplots(figsize=(10,6), dpi=500)
ax1.set_xlabel('Time (sec)', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm05, label='Change in binding motor')
ax1.legend(loc='upper left')
ax1.set_ylim(bottom=0)
ax2 = ax1.twinx()
ax2.set_xlabel('Time (sec)', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(nbm05.mean(), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')
ax2.set_ylim(bottom=0)
plt.title('Binding motors [R = 0.5]', fontdict=font)
plt.grid()
plt.show()
fig.savefig('actin_binding_grinding_r0.5.svg',bbox_inches='tight', format='svg')
#print(colored('Ave_vel:', 'yellow', attrs=['reverse', 'blink']), \
# colored(Av_vel, 'yellow', attrs=['reverse', 'blink']))
#print(colored('VSD:', 'yellow', attrs=['reverse', 'blink']), \
# colored(vSD, 'yellow', attrs=['reverse', 'blink']))
print(colored('Mean velocity, time[0:0.3] = ', 'yellow', attrs=['reverse', 'blink']), \
colored(nbm05.mean(), 'yellow', attrs=['reverse', 'blink']))
print(colored('bVSD:', 'yellow', attrs=['reverse', 'blink']), \
colored(nvSD, 'yellow', attrs=['reverse', 'blink']))
nbm05 = pd.concat([tym, bm05_, v1], axis=1)
nbm05.columns=['time','b_m','vel']
print(nbm05)
Dx_tip = np.diff(bmd06['x']); Dy_tip = np.diff(bmd06['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(10*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.around(np.linspace(0,3,30,endpoint=True), decimals=1)
fig, ax1 = plt.subplots(figsize=(10,6), dpi=500)
ax1.set_xlabel('Time (sec)', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm06, label='Change in binding motor')
ax1.legend(loc='upper left')
ax1.set_ylim(bottom=0)
ax2 = ax1.twinx()
ax2.set_xlabel('Time (sec)', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')
ax2.set_ylim(bottom=0)
plt.title('Binding motors [R = 0.6]', fontdict=font)
plt.grid()
v1 = v.reshape(30,1)
v1 = pd.DataFrame(v1)
tym = pd.DataFrame(time)
bm06_ = bm06.reset_index(drop=True)
nbm06 = pd.concat([tym, bm06_, v1], axis=1)
nbm06.columns=['time','b_m','vel']
#print(nbm06)
plt.show()
fig.savefig('actin_binding_grinding_r0.6.svg',bbox_inches='tight', format='svg')
print(colored('Ave_vel:', 'yellow', attrs=['reverse', 'blink']), \
colored(Av_vel, 'yellow', attrs=['reverse', 'blink']))
print(colored('VSD:', 'yellow', attrs=['reverse', 'blink']), \
colored(vSD, 'yellow', attrs=['reverse', 'blink']))
nbm06 = pd.concat([tym, bm06_, v1], axis=1)
nbm06.columns=['time','b_m','vel']
print(nbm06)
Dx_tip = np.diff(bmd07['x']); Dy_tip = np.diff(bmd07['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(10*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.around(np.linspace(0,3,30,endpoint=True), decimals=1)
fig, ax1 = plt.subplots(figsize=(10,6), dpi=500)
ax1.set_xlabel('Time (sec)', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm07, label='Change in binding motor')
ax1.legend(loc='upper left')
ax1.set_ylim(bottom=0)
ax2 = ax1.twinx()
ax2.set_xlabel('Time (sec)', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')
ax2.set_ylim(bottom=0)
plt.title('Binding motors [R = 0.7]', fontdict=font)
plt.grid()
v1 = v.reshape(30,1)
v1 = pd.DataFrame(v1)
tym = pd.DataFrame(time)
bm07_ = bm07.reset_index(drop=True)
nbm07 = pd.concat([tym, bm07_, v1], axis=1)
nbm07.columns=['time','b_m','vel']
#print(nbm07)
plt.show()
fig.savefig('actin_binding_grinding_r0.7.svg',bbox_inches='tight', format='svg')
print(colored('Ave_vel:', 'yellow', attrs=['reverse', 'blink']), \
colored(Av_vel, 'yellow', attrs=['reverse', 'blink']))
print(colored('VSD:', 'yellow', attrs=['reverse', 'blink']), \
colored(vSD, 'yellow', attrs=['reverse', 'blink']))
nbm07 = pd.concat([tym, bm07_, v1], axis=1)
nbm07.columns=['time','b_m','vel']
print(nbm07)
Dx_tip = np.diff(bmd08['x']); Dy_tip = np.diff(bmd08['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(10*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.around(np.linspace(0,3,30,endpoint=True), decimals=1)
fig, ax1 = plt.subplots(figsize=(10,6), dpi=500)
ax1.set_xlabel('Time (sec)', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm08, label='Change in binding motor')
ax1.legend(loc='upper left')
ax1.set_ylim(bottom=0)
ax2 = ax1.twinx()
ax2.set_xlabel('Time (sec)', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')
ax2.set_ylim(bottom=0)
plt.title('Binding motors [R = 0.8]', fontdict=font)
plt.grid()
v1 = v.reshape(30,1)
v1 = pd.DataFrame(v1)
tym = pd.DataFrame(time)
bm08_ = bm08.reset_index(drop=True)
nbm08 = pd.concat([tym, bm08_, v1], axis=1)
nbm08.columns=['time','b_m','vel']
#print(nbm08)
plt.show()
fig.savefig('actin_binding_grinding_r0.8.svg',bbox_inches='tight', format='svg')
print(colored('Ave_vel:', 'yellow', attrs=['reverse', 'blink']), \
colored(Av_vel, 'yellow', attrs=['reverse', 'blink']))
print(colored('VSD:', 'yellow', attrs=['reverse', 'blink']), \
colored(vSD, 'yellow', attrs=['reverse', 'blink']))
nbm08 = pd.concat([tym, bm08_, v1], axis=1)
nbm08.columns=['time','b_m','vel']
print(nbm08)
Dx_tip = np.diff(bmd09['x']); Dy_tip = np.diff(bmd09['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(10*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.around(np.linspace(0,3,30,endpoint=True), decimals=1)
fig, ax1 = plt.subplots(figsize=(10,6), dpi=500)
ax1.set_xlabel('Time (sec)', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm09, label='Change in binding motor')
ax1.legend(loc='upper left')
ax1.set_ylim(bottom=0)
ax2 = ax1.twinx()
ax2.set_xlabel('Time (sec)', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')
ax2.set_ylim(bottom=0)
plt.title('Binding motors [R = 0.9]', fontdict=font)
plt.grid()
v1 = v.reshape(30,1)
v1 = pd.DataFrame(v1)
tym = pd.DataFrame(time)
bm09_ = bm09.reset_index(drop=True)
nbm09 = pd.concat([tym, bm09_, v1], axis=1)
nbm09.columns=['time','b_m','vel']
#print(nbm09)
plt.show()
fig.savefig('actin_binding_grinding_r0.9.svg',bbox_inches='tight', format='svg')
print(colored('Ave_vel:', 'yellow', attrs=['reverse', 'blink']), \
colored(Av_vel, 'yellow', attrs=['reverse', 'blink']))
print(colored('VSD:', 'yellow', attrs=['reverse', 'blink']), \
colored(vSD, 'yellow', attrs=['reverse', 'blink']))
nbm09 = pd.concat([tym, bm09_, v1], axis=1)
nbm09.columns=['time','b_m','vel']
print(nbm09)
Dx_tip = np.diff(bmd10['x']); Dy_tip = np.diff(bmd10['y'])
DD=np.sqrt((Dx_tip**2)+(Dy_tip**2))
v=DD/(10*dt); Av_vel = np.mean(v)
vSD=np.sum(((v-Av_vel)**2)/(np.size(v)-1)); vSD=np.sqrt(vSD)
time = np.around(np.linspace(0,3,30,endpoint=True), decimals=1)
fig, ax1 = plt.subplots(figsize=(10,6), dpi=500)
ax1.set_xlabel('Time (sec)', fontdict=font)
ax1.set_ylabel('Number of binding motors', fontdict=font)
ax1.plot(time,bm10, label='Change in binding motor')
ax1.legend(loc='upper left')
ax1.set_ylim(bottom=0)
ax2 = ax1.twinx()
ax2.set_xlabel('Time (sec)', fontdict=font)
ax2.set_ylabel('Velocity ($\mu m/sec$)', fontdict=font)
ax2.plot(time,v, 'g', label='Change in velocity')
ax2.axhline(np.mean(v), linestyle='--', color='g', label='Mean velocity')
ax2.legend(loc='upper right')
ax2.tick_params('y', colors='g')
ax2.set_ylim(bottom=0)
plt.title('Binding motors [R = 1.0]', fontdict=font)
plt.grid()
v1 = v.reshape(30,1)
v1 = pd.DataFrame(v1)
tym = pd.DataFrame(time)
bm10_ = bm10.reset_index(drop=True)
nbm10 = pd.concat([tym, bm10_, v1], axis=1)
nbm10.columns=['time','b_m','vel']
#print(nbm10)
plt.show()
fig.savefig('actin_binding_grinding_r1.0.svg',bbox_inches='tight', format='svg')
print(colored('Ave_vel:', 'yellow', attrs=['reverse', 'blink']), \
colored(Av_vel, 'yellow', attrs=['reverse', 'blink']))
print(colored('VSD:', 'yellow', attrs=['reverse', 'blink']), \
colored(vSD, 'yellow', attrs=['reverse', 'blink']))
nbm10 = pd.concat([tym, bm10_, v1], axis=1)
nbm10.columns=['time','b_m','vel']
print(nbm10)
columns = ['r', 'b', 'atp', 'md', 'v_av', 'v_sd']
df = pd.read_csv('bm_grinding_vel.csv', names=columns)
x=df['r']; y=df['v_av']; dev=df['v_sd']
plt.figure(figsize=(10,6), dpi=500)
plt.errorbar(x,y,yerr=dev, ecolor='b', capsize=7, \
mec='blue', color='green', marker='D', linewidth=2, \
markersize=7, label='Velocity deviation')
plt.xlabel('Specie ratio ($\%$)', fontdict=font); plt.ylabel('Velocity ($\mu m/sec$)', fontdict=font)
plt.title('Actin Filament Grinding Movement [ATP = 2000, MD = 3000]'); plt.legend(loc='upper left'); plt.grid()
plt.savefig('actin_bm_grinding.svg',bbox_inches='tight', format='svg')
print(colored('Data for R, Velocity and Deviation:', 'yellow', attrs=['reverse', 'blink']))
print(df)